Eigenstate Structure in Graphs and Disordered Lattices 

L. Kaplan* 

Institute for Nuclear Theory and Department of Physics 
University of Washington, Seattle, Washington 98195 



o 
o 

(N 

c 

(N 



U 
d 



> 
oo 

O 

o 
o 

• 1— I 

> 
X 



We study wave function structure for quantum graphs in 
the chaotic and disordered regime, using measures such as 
the wave function intensity distribution and the inverse par- 
ticipation ratio. The result is much less ergodicity than ex- 
pected from random matrix theory, even though the spectral 
statistics are in agreement with random matrix predictions. 
Instead, analytical calculations based on short-time semiclas- 
sical behavior correctly describe the eigenstate structure. 

Quantum graphs, known also as network models, have 
been used successfully for many years as simple dynam- 
ical systems in which to study complex wave behavior. 
For appropriate parameter values, graphs can be made 
to display generic chaotic, disordered, or integrable mo- 
tion, and at the same time the quantum mechanics of 
these systems has the simplifying advantage of being 
semiclassically exact. Originally, graphs were developed 
as models for electrons moving between atoms in an or- 
ganic molecule 0]]. Recently, graph models have been 
used to study issues as diverse as Anderson localization 
within the context of periodic orbit theory || , the spatial 
distribution and transport properties of persistent cur- 
rents H|, Aharonov-Bohm conductance modulations in 
GaAs/GaAlAs networks spectral statistics and the 
trace formula in chaotic systems spectral determi- 
nants, with applications to thermodynamic and trans- 
port properties of mesoscopic networks [||, and chaotic 
scattering and resonance behavior J?J. A discussion of 
the earlier history of quantum graphs can be found in 
Ref. [jjj, which also provides an extensive discussion of 
the model. 

Though substantial work now exists on the spectral 
and scattering properties of quantum graphs, and also on 
their large-scale localization behavior, surprisingly little 
attention has been paid so far to the detailed wave func- 
tion structure of this paradigmatic quantum system. In 
this paper we begin to address questions relating to the 
statistics of wave functions on graphs and their relation to 
the underlying classical structures in the system. In the 
process, we examine the relationship between short-time 
and long-time effects on stationary behavior, making con- 
nections to recent work on other simple quantum chaotic 
problems, including quantum maps, Sinai billiards, Buni- 
movich stadia, tunneling in double wells, conductance 
through chaotic quantum dots, and many-body systems 
with random two-body interactions Q. 



A quantum graph consists of V vertices connected by 
B bonds, each of which is modeled as a one-dimensional 
wire of length Lj (j = 1 . . . B) along which the wave 
propagates freely with zero potential. At each of the 
i = 1 . . . V vertices, Vi > 2 bonds will meet (note that 
JT Vi = 2B), and one must impose wave function conti- 
nuity and the current conservation condition 
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where the sum is over all Vi bonds j that meet at ver- 
tex i and ^ij(x) is the wave function in the bond j, 
with x — corresponding to the vertex i. A; is a free 
parameter associated with the height of the effective po- 
tential at vertex i, and allows one to interpolate between 
Neumann (A.; = 0) and Dirichlet (A^ — > oo) boundary 
conditions. One sees easily that the form of Eq. (Q) is 
the only one consistent with an s-wave pointlike scat- 
terer at the vertex; it can be obtained formally using 
self-adjoint extension theory. Because of time reversal 
invariance, the wave function in each bond can be writ- 
ten as ^ij(x) — a^e lkx + af^j* : e ~ lkx , for an eigenstate 
at energy k 2 /2m. Thus, the distribution of wave function 
intensities in the system can be completely characterized, 
in the limit kL — > oo, by the distribution of the quantities 
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In this limit, the simple normalization condition 
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is set to 
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ensures that the mean bond intensity 
unity. 

We begin our analysis with a simple one-dimensional 
"ring graph" Q where the V vertices are arranged in a 
circle and each vertex i is connected by bonds to neigh- 
boring vertices i — v/2 . . . i — l,i + 1 . . .i + v/2. v, the 
number of bonds meeting at every vertex, is known as the 
valency (here taken to be constant over the entire graph) . 
Setting all the vertex potentials Xi — leads to the max- 
imum possible delocalization in the graph (the opposite 
limit, Xi — > oo would instead produce eigenstates local- 
ized on individual bonds), while non-integrability is en- 
sured through randomness in the bond lengths Lj. In 
our calculation we take the bond lengths to be uniformly 
distributed in an interval [1— SL, 1+5L); because the scat- 
tering matrices depend only on kLj mod 2tt, all choices 
of SL arc equivalent as long as SL ^ fc -1 . The eigen- 
states of the system may be obtained by examining a 
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V x V secular matrix h(k) [||; k is an eigenvalue when- 
ever det h(k) = 0, and the associated null vector cor- 
responds to an eigenvector of the graph. In fact, small 
but non-zero singular values e of h(k) can easily be seen 
to correspond to eigenvectors of the same system with 
slightly perturbed potentials Aj — * Aj — fee, so sufficiently 
small singular values (|e| < fc _1 A) can also be used to 
produce eigenstates. This method allows the collection 
of several independent wave functions at a given value 
of k for any given realization of the disorder ensemble. 
One easily checks that collecting one or more wave func- 
tions at a given value of k has no discernible effect on the 
resulting wave function statistics. 

Comparing the Heisenberg time (h divided by the 
mean level spacing, at which individual levels are re- 
solved) with the time required classically to diffuse over 
the entire system, or alternatively making an analogy 
with banded random-matrix behavior, we see that the 
condition for avoiding localization in our system is v 2 3> 
V. We note that the localization condition is fc— and 
^—independent and depends only on the classical graph 
geometry. Increasing the valency v for a fixed system size 
V, one easily observes a transition from localized to de- 
localized behavior, which can be detected either by look- 
ing at the change in level spacing statistics (from Poisson 
to GOE) or at the change in the wave function intensity 
correlation (from strongly negative to near zero) between 
distant points on the graph. Either method confirms the 
expected scaling behavior for the transition. One can 
therefore take the large- volume limit V — > oo, where sta- 
tistical behavior is expected, while easily satisfying the 
derealization condition y/V < v < V — 1. 

In the delocalized regime, full information about wave 
function statistical behavior is contained in the distribu- 
tion of bond intensities \a,j | 2 and their correlations (note 
that we may freely drop the i index as it is immaterial 
which of the two endpoints we take to be the beginning 
of bond j). It is convenient to introduce a simple one- 
number measure of wave function ergodicity, the inverse 
participation ratio (IPR): 



I=(\a[ k) 
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where the averaging is performed over the B — Vv/2 
bonds j and over a disorder ensemble, at a fixed value of 
k. Of course, averaging over nearby values of k may also 
be done. It is often useful to introduce a local version of 
the IPR, 2j, where the bond j is fixed; then 1 = (2j), 
averaging over all bonds. The IPR is the first nontrivial 
moment of the intensity distribution P(|ay | 2 ) (we recall 
that the mean intensity has been normalized to unity), 
and can range from 1 in the maximally ergodic case where 
all intensities are equal up to a maximum value of B in 
the case of perfect wave function localization on individ- 
ual bonds. Random matrix theory or the random vector 
hypothesis would predict Gaussian random fluctuations 
in the complex coefficients a^, thus 1=2. An enhance- 
ment of the IPR above this baseline value indicates a 



deviation from ergodicity in the local wave function be- 
havior. 

The key theoretical idea discussed and applied in sev- 
eral recent works is that wave function intensities in a 
complex system can often be conveniently separated into 
a product of short-time and long-time parts: 
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Pj{E{k)) x r jk 
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Here Pj(E) is a smooth local density of states (known 
alternatively as the strength function) on the bond j at 
energy E, obtained as the Fourier transform of the short- 
time part of the autocorrelation function 
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while rjk is obtained (formally) by Fourier transforming 
the long-time behavior, at times t ~ T cuto ff and larger. 
The decomposition is useful because in many situations 
the short-time return amplitude A:j hort (i) has a known 
approximate analytical expression, which can be trans- 
formed to obtain pj(E). On the other hand the long- 
time return amplitude in a chaotic or disordered system 
is given by a convolution of the short-time behavior with 
a sum of exponentially many contributions, and thus 
may be regarded as random variable: 



{r jk ) = 1 
(r jk r 3 >k>) =l + (F - l)5jf5 k k' 
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where the statistical average is performed over an ap- 
propriate ensemble. (If rjk is the square of a complex 
Gaussian random variable, then F = 2.) Because the 
smooth local spectral density pj (E) and the fluctuations 
rjk are associated with distinct time scales (before and 
after the mixing time, respectively, in a chaotic system), 
the two quantities are regarded as statistically indepen- 
dent. Thus, for example, the local IPR can be written 



(P 2 j)(r%) = (P 2 3 ) 



F. 



(7) 



where (p 2 ), the second moment of the smooth local den- 
sity function, is proportional to the sum of short-time 
return probabilities |Aj hort (i)| 2 . This formalism has suc- 
cessfully been used to quantitatively study scars of un- 
stable periodic orbits and related phenomena in billiards, 
in smooth potential wells, and in many-body interacting 
systems. 

To apply these ideas to the ring graphs, we focus on 
one (arbitrary) bond j connecting vertices 1,2. An initial 
wave packet launched in this bond moving from 1 towards 
2 will have a probability 
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of being transmitted into one of the other v — 1 bonds 
meeting at vertex 2; the remaining part of the wave 
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packet then gets reflected back into the original bond || . 
To begin with, we set A = at all vertices for sim- 
plicity, and find that the reflected probability is P rc fl = 
1 — (v — l)Ptrans = 1 — 4(v — 1) / ' v 2 . The process is re- 
peated at vertex 1, and the remaining probability P r g fl 
travels again the path taken by the original wave packet, 
leading to a nontrivial contribution to the return prob- 
ability. We may iterate this process until almost all of 
the initial probability to be in the bond j has decayed 
(after 0(1/ v) bounces for v 3> 1), and find that the sum 
of return probabilities behaves as 



dt 
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(note that we always sum return probabilities over both 
positive and negative times). To leading order in v, 
we therefore obtain (p 2 ) = v/4 — 0(1) for the short- 
time factor. Taking into account intermediate-time re- 
currences (where the wave is transmitted at either vertex 
into an adjoining bond and is subsequently transmitted 
back into the bond j) cancels the 0(1) term, leading to 
(Pj ) = v /4 + Oiv- 1 ), and thus 
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FIG. 1. Observed IPR vs. that predicted by the short-time 
theory, Eq. (0), for ring graphs with parameter values 
7 < v < 15 and 15 < V < 33. Data points represented 
by + signs for v — 7, crosses for v = 9, stars for v = 11, 
empty squares for v = 13 and full squares for v = 15. 

We note that in Eq. ([n]) we have included the leading 
0(h) correction to the semiclassical answer (for graph 
systems V -1 plays the role of an effective h, and b is an 
undetermined constant). The long-time factor F can be 
obtained directly by measuring the mean square value of 
\aj\ 2 /pj(E(k)), i.e. of the bond intensity normalized by 
the analytically computed short-time spectral envelope. 
For a wide range of values for the system size V and the 
valency v, the result 4 < F < 5 is obtained, supporting 
the conjecture of independence of short-time and long- 
time fluctuations in the spectrum (we note F > 2 means 



that the long-time fluctuations are super-Gaussian). Nu- 
merical data can then be used to fit the coefficient b of the 
subleading semiclassical correction; the fit is quite good 
as can be seen in Fig. [j]. As expected, the IPR is pri- 
marily a function of valency v, with volume-dependence 
being a higher-order effect as long as the derealization 
condition v 2 > V is satisfied. The wave function statis- 
tics clearly deviate strongly from random matrix expecta- 
tions as v becomes large. This is despite the fact that the 
level spacing statistics of this system are well-predicted 
by random matrix theory, indicating an absence of strong 
localization. 

In a one-dimensional system it is of course impossible 
to take the semiclassical (large-volume) limit for fixed v 
while staying in the delocalized regime. It is therefore of 
interest to consider higher-dimensional systems, such as a 
d-dimensional cubic lattice, with v = 2d. In the absence 
of vertex potentials (A = 0) the above analysis still ap- 
plies. What happens when we introduce disorder into the 
system via the potentials Xi in addition to the disorder 
already present in the bond lengths? Let the A^ be inde- 
pendent and distributed for large A^ in accordance with a 
power law T{\i) ~ A^ Q , for A < Xi < oo (with a > 1). 
We now claim that the tail of the IPR distribution will be 
strongly modified by the rare events where a strong po- 
tential A is present on both sides of a given bond j with 
endpoints 1 and 2. Indeed, we easily see that Eq. (j|) 
for transmission probability reduces to Ptrans ~ A~ 2 for 
strong A. Clearly the weaker of the two potentials Ai 
and A2 will dominate the escape rate. The short-time 
enhancement factor for the local IPR is proportional to 
the inverse of the escape rate, i.e. Tj ~ min(A^, X 2 ), and 
thus we have the prediction 
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for 1 <§; Tj <C V, modifying the exponential falloff pre- 



dicted by random matrix theory, 
the intensity distribution becomes 



V(\a\ 2 )~X 
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Similarly, the tail of 
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for 1 < 



<^ V, in contrast with the exponential 



Porter-Thomas prediction valid for a system satisfying 
random matrix statistics. 
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FIG. 2. The tail of the wave function intensity distribution 
is plotted for a three-dimensional lattice with random poten- 
tials A (see text). From top to bottom, the three solid data 
curves are for exponent a = 1.2, 1.5, and 1.8. The corre- 
sponding theoretical power laws from Eq. ( |l2j ) are plotted as 
dashed lines, with the RMT prediction appearing as a dotted 
curve for comparison. 

The result of Eq. ( jl2| ) is confirmed in Fig. ||, where 
we have used an ensemble of 37 3 lattices with one out of 
every D = 11 sites randomly chosen to contain a scat- 
terer (the particle executes free motion while traveling 
between the scatterer sites, so that the classical motion 
is diffusive with mean free path D). The total number of 
scatterers is then V = 4604, and we indeed see that the 
data curves tend quickly to zero for \a\ 2 > V. We have 
checked that with increasing number of scatterers V at 
fixed a, the intensity distribution curves keep following a 
given power law behavior for larger and larger intensity, 
before eventually dropping off to zero. Similarly, we have 
checked that varying the diffusion constant I? by a factor 
of 2 does not significantly affect the intensity distribution 
for fixed V, as would have been expected for a diffusion- 
dominated log-normal tail, logV ~ — Z?log 2 (|a| 2 ) [[To|| . 
We see instead that the tail of the intensity distribution 
is dominated entirely by the short-time system behavior. 

Finally, we return to the volume-averaged IPR X, the 
simplest overall measure of the degree of wave function 
localization. From Eq. (p~T|) , we easily obtain in the large- 
volume limit 



X — C a A 



(13) 



for 1 < a < 2. As we can see in Table I, this result 
compares favorably to the numerically computed value 
over a range of volumes V and for different classical mean 
free paths D and minimum potentials Ao- None of the 
IPR's are close to the random matrix prediction 2 = 2. 
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TABLE I. Numerically obtained inverse participation ratio 
T for a disordered three-dimensional lattice with exponent 
a = 1.5, compared with the prediction of Eq. ( |l3| ) for various 
numbers of scatterers V, mean free paths D, and minimum 
potential values Ao- The constant used is C a = 3.0. 



In conclusion, we have seen that various aspects of 
wave function structure in chaotic and disordered quan- 
tum graphs are adequately described using the analyti- 
cally known short-time behavior of the system, specifi- 
cally the return amplitude at short times and its Fourier 
transform. On the other hand, random matrix theory 
completely fails to describe the wave function structure, 
even though its predictions are good for the spectral 
statistics. This failure can be understood as resulting 
from the omnipresence of short periodic orbits in graphs, 
in contrast with the situation prevailing in most other 
chaotic systems. 

Several very valuable discussions with T. Kottos are 
gratefully acknowledged. This work was supported by 
the DOE under grant No. DE-FG03-00-ER41132. 



[1 
[2 
[3. 

[4: 
[5] 

[6] 

[7 

[8] 
[9] 



M. J. Richardson and N. L. Balazs, Ann. Phys. (N. Y.) 
73, 308 (1972). 

H. Schanz and U. Smilansky, Phys. Rev. Lett. 84, 1427 
(2000). 

M. Pascaud and G. Montambaux, Phys. Rev. Lett. 82, 
4512 (1999). 

J. Vidal, G. Montambaux, and B. Doucot, Phys. Rev. E 
62, R16294 (2000). 



H. Schanz and U. Smilansky, |arXiv:chao-dyn/9904007 . 

E. Akkermans, A. Comtet, J. Desbois, G. Montambaux, 

and C. Texier, J. Phys. A 33, L63 (2000). 

T. Kottos and U. Smilansky, Phys. Rev. Lett. 85, 968 

(2000). 

T. Kottos and U. Smilansky, Ann. Phys. 274, 76 (1999). 
L. Kaplan and T. Papenbrock, Phys. Rev. Lett. 84, 4553 
(2000); L. Kaplan and E. J. Heller, Phys. Rev. E 62, 409 
(2000); L. Kaplan, Phys. Rev. E 62, 3476 (2000); W. 
E. Bies, L. Kaplan, M . R. Haggerty, and E. J. Heller, 
arXiv:nlin.CD/0004024j; W. E. B ies, L. Kaplan, and E. 
J. Heller, |axXiv:nlm.CD/0007037j ; L. Kaplan, Phys. Rev. 



Lett. 80, 2582 (1998). 
[10] V. I Fal'ko and K. B. Efetov, Phys. Rev. B 52, 17413 
(1995); I. E. Smolyarenko and B. L. Altshuler, Phys. Rev. 
B 55, 10451 (1997). 



4 



